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Abstract 



■ ■ We investigate the question of possible integrability of classical string motion in curved p- 

, brane backgrounds. For example, the D3-brane metric interpolates between the flat and the 

AdSc, X regions in which string propagation is integrable. We find that while the point- 
like string (geodesic) equations are integrable, the equations describing an extended string in 
' the complete D3-brane geometry are not. The same conclusion is reached for similar brane 

■ intersection backgrounds interpolating between flat space and AdSk x S**-'. We consider, in 

, particular, the case of the NS 5-brane - fundamental string background. To demonstrate non- 

I integrability we make a special "pulsating string" ansatz for which the string equations reduce 

to an effective one-dimensional system. Expanding near this simple solution leads to a linear 
differential equation for small fluctuations that cannot be solved in quadratures, implying non- 
integr ability of the original set of string equations. 
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1 Introduction 



Strings in curved space are described by 2d sigma models with complicated non-linear equations of 
motion. This precludes, in most detailed understanding of their dynamics. Integrability 

selects a subclass of models or target space backgrounds for which one may hope to come as close 
as possible to the situation in flat space, i.e. to have a complete quantitative description of classical 
string motions and then also of the quantum string spectrum. 

A prominent maximally symmetric 10-dimensional example is provided by the superstring in the 
AdS^ X space-time: much progress towards its solution, that was achieved recently [1], is based 
on its integrability. One may wonder if integrability can also apply to string motion in closely 
related but less symmetric p-brane backgrounds [2-4], e.g. the D3-brane background [4,5] which is 
a 1-parameter (D3-brane charge Q) interpolation between flat space (Q = 0) and the AdS^ x 
space {Q = oo). While string theory is integrable in these two limiting cases, what happens at finite 
Q is a priori an open question. This is the question we are going to address in this paper. 

Since quantum integrability of a bosonic model or integrability of its superstring counterpart is 
essentially implied (leaving aside some exceptional cases or anomalies) by the integrability of the 
corresponding classical bosonic sigma model the study of the latter is the first step. Following the 
same method as used previously in [6,7] we will demonstrate non-integrability of classical extended 
string motion in the D3-brane background (particle, i.e. geodesic, motion is still integrable). We 
will also reach similar conclusions for other p-brane backgrounds that interpolate between fiat space 
and integrable AdSn x backgrounds. 

As there is no general classification of integrable 2d sigma models, let us start with a brief 
survey of some known classically integrable bosonic models describing string propagation in curved 
backgrounds. A major class of such models have a target space metric Gmn of a symmetric space, 
including spheres [8] (and other related constant curvature spaces [9]), group spaces (i.e. principal 
chiral model) [10] and various other G/H coset models [10-12]. Integrability is also preserved by 
some anisotropic deformations of the group space [13,14], the 2-sphere [15] and the 3-sphere [16-18] 
models. Attempts to find "non-diagonal" generalisations of the principal chiral model by studying 
conditions for the existence of the corresponding zero curvature (Lax) representation were made 
in [19,20] and also in [21]. 

Allowing for the parity-odd antisymmetric tensor {Bmn) coupling leads, in particular, to the 
WZW and related gauged WZW models [22, 23]. The generalised principal chiral model with 
the WZ term having an arbitrary coefficient is also integrable [24-26]. Some other "anisotropic" 
integrable models with 3d target space were constructed in [18] and new examples of integrable 
coset-type models with extra WZ coupling were found in [26]. 

There are also examples of integrable pp-wave models related to (massive) integrable 2d models 
in light-cone gauge, see e.g. [27-31]. 
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Given an integrable model one can generate a family of other integrable models (i.e. find new 
"integrable" target space backgrounds) by performing transformations that preserve the classical 2d 
equations of motion, e.g. using 2d abelian duality (T-duality) or non-abelian duality, combined with 
field redefinitions (coordinate transformations). Examples include, in particular, various gaugings 
or marginal deformations of WZW and related models, see e.g. [32-36]; for a more recent example 
related to AdS^ x see [37].^ For discussions of an interplay between T-duality and integrability 
see [21,40]. 

Given a particular string sigma model, to prove its integrability one is to find a zero curvature 
representation implying its equations of motion. There is no general method for accomplishing this 
and the full set of necessary conditions for the existence of a Lax pair is not known. For example, 
the presence of a non-abelian isometry group is not required.^ Among the necessary conditions 
should, of course, be that any consistent truncation of the 2d equations of motion to a Id set of 
equations has to represent an integrable mechanical system. For example, rigid string motion on a 
sphere or AdS space are described by an integrable Neumann- Rosochatius model [41]. 

In particular, the point-like string, i.e. geodesic, motion should be integrable. The geodesic 
motion in a general class of higher-dimensional (rotating) black hole backgrounds is known to be 
integrable [42].^ It is thus not surprising that geodesies are also integrable in p-brane backgrounds 
that are discussed below. 

To disprove integrability of a given sigma model it is therefore sufficient to show that there exists 
at least one Id truncation of its equations of motion that is not an integrable system of ordinary 
second-order differential equations. This can be done by demonstrating that the corresponding 
motion is chaotic, i.e. by showing that the equations for small variations around solutions in phase 
space cannot be solved in quadratures [48]. We shall review this method in section 2 (and Appendix 
A) and illustrate it on a simple example of the non-integrable two-particle system with the potential 
V{xi,X2) = \xlx'l. 

This approach was recently used to show that string motion in AdS^ x T^'"?, AdS^ x Y'P''^ [6,7] and 
"confining" supergravity backgrounds [49] is not integrable.^ Here we shall follow a similar method 

^One may also find new models by taking orbifolds of known integrable models (see, e.g., [38]), though division 
over discrete subgroups may also lead to a breakdown of integrability (of. [39]). 

■^Indeed, generic sigma models obtained from gauged WZW models by solving for the 2d gauge field do not have 
isometries (or have only abelian isometries) and yet, following from a gauged WZW action, they must be integrable. 

^Also, a particular stationary string motion in AdS-Kerr-NUT backgrounds is integrable as well [43]. 

*In addition, if the general string motion is integrable, then motion of a subclass of strings with large momentum 
along a null geodesic or, equivalently, string motion in a pp-wave background which is a Penrose limit of a given 
curved background should of course be integrable as well. For the p-brane backgrounds that we are going to consider 
string motion in the pp-wave backgrounds representing Penrose limits turn out to be solvable (see e.g. ]44-47]) as 
in the light-cone gauge description the string moves in a quadratic potential. However, since the coefficient of this 
potential happens to depend on r in general, the corresponding models do not appear to be integrable. 

^Spaces like T^'^ and Y^''^ are special in that they have isometries, preserve supersymmetry and have integrable 
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to study integrability of strings in curved backgrounds that interpolate between two integrable limits 
- flat space and AdSn x S"" x T''. 

The main cases include the D3-brane, the D5-D1 background and the four D3-brane intersection 
that interpolate between flat space and AdS^ x , AdS^ x x and AdS2 x 5^ x respectively. 
While geodesic motion in these backgrounds is integrable, making a special ansatz for string motion 
we will find that the resulting Id system of equations is not integrable by applying the variational 
non-integrability technique for Hamiltonian systems as used in [7] (section 3). 

In section 4 we shall study the NS5-F1 background and show that integrability is absent for 
generic values of the two charges Qi and Q^, but it is restored in the limit when — ?• oo. 

In Appendix A we also give a brief summary of the ideas from differential Galois theory on 
which the non-integrability techniques for Hamiltonian systems are based on. In Appendix B we 
shall demonstrate non-integrability of string propagation in the pp-wave geometry T-dual to the 
fundamental string background [30,51]. 

2 On the non-integrability of classical Hamiltonian systems 

Integrability of classical Hamiltonian systems is closely related to the behaviour of variations around 
phase space curves. Based on this observation, ref. [52] obtained necessary conditions for the 
existence of additional functionally independent integrals of motion. These conditions are given 
in terms of the monodromy group properties of the equations for small variations around phase 
space curves. Using differential Galois theory, refs. [48,53,54] improved these results by showing 
that integrability implies that the identity component of the differential Galois group of variational 
equations normal to an integrable plane of solutions must be abelian (see Appendix A for some 
details) . 

This necessary integrability condition allows one to show non-integrability from the properties of 
the Galois group. However, determining the Galois group can be difficult and therefore one usually 
takes a slightly different route in analysing the normal variational equation (NVE). 

One considers a special class of solutions of the NVE which are functions of exponentials, loga- 
rithms, algebraic expressions of the independent variables and their integrals and which are known 
as Liouvillian solutions [55]. The existence of such solutions is equivalent to the condition that the 
identity component of the Galois group is solvable (see, for example, theorem 25 in [56]). This 
in turn implies that if the NVE does not admit Liouvillian solutions then is non-solvable (and 
thus non-abelian) leading to the conclusion of non-integrability. 

For Hamiltonian systems the NVE is a second order linear homogeneous differential equation and 

geodesies (in fact, super-integrable - with more conserved quantities than degrees of freedom [50]). Yet, the corre- 
sponding 2d sigma models fail to be integrable. 



3 



for such equations with rational coefficients Liouvillian solutions can be determined by the Kovacic 
algorithm [57] which fails if and only if no such solutions exist. Therefore, whenever no Liouvillian 
solutions of the NVE are found by the Kovacic algorithm one can deduce non-integrability of a 
Hamiltonian system. 

An important subtlety here is that we first need to algebrize the NVE, i.e. rewrite it as a 
differential equation with rational coefficients. This can be done by means of a Hamiltonian change 
of the independent variable and it was shown in [58] that the identity component of the Galois 
group is preserved under this procedure. 

In summary, if the algebrized NVE is not solvable in terms of Liouvillian functions or, equivalently, 
if the component of the quadratic fiuctuation operator normal to an integrable subsystem does not 
admit Liouvillian zero modes, the original system is non-integrable. This is consistent with the usual 
definition of integrability in the sense of Liouville which, for Hamiltonian systems, implies that the 
equations of motion should be solvable in quadratures, i.e. in terms of Liouvillian functions. 

The main steps to prove non-integrability of a Hamiltonian system are thus: 

• Choose an invariant plane of solutions. 

• Obtain the NVEs, i.e. the variational equations normal to the invariant plane. 

• Check that the algebrized NVEs have no Liouvillian solutions using the Kovacic algorithm. 
This approach was recently used in [7,49] to show the non-integrability of string motion in several 

curved backgrounds and we shall also follow it here. 



Let us first explain what the algebraization procedure of the last step means. A generic NVE is 
a second order differential equation of the form 

fi + q{t)iQ + r{t)ri = . (2.1) 

A change of the variable t — t- x{t) is called Hamiltonian if x{t) is a solution of the Hamiltonian 

2 

system H{p, x) = ^ + V{x) where V{x) is some potential. Then x satisfies the first integral 
equation -|- V{x) = h, implying that x (and also x) is a function of x only. Changing the 
variable t — t- x(t) in (2.1) we obtain (here ' = ^) 

„ /X q(t(x))\ , r(t(x)) 
r]" + +^\^V = . 2.2 

Vx^ X / x^ 

We should now require that the coefficients in this equation are rational functions of x. 

Let us consider the following simple example of a Hamiltonian system with the canonical variables 
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{xi,pi) (i = l,2) [59] 



H = 1{pI+pI) + 1x1x1 (2.3) 



A choice of an invariant plane, referred to above, is P = {{xi,X2,pi,P2) ■ X2 = P2 = 0} along which 
the solution curves are xi(t) = kt + const, pi = k, k = \/2H. The system along the plane P with 
the coordinates {xi,pi) is integrable since the only required constant of motion is provided by the 
Hamiltonian. 

To determine the direction normal to the solution curves of this integrable subsystem we can 
use the Hamiltonian vector field along the curves X^^ = k, X^^ = 0, X^^ = 0, = 0. Thus 

the normal direction is along X2 and p2 and expanding around X2 = with 6x2 = rj we obtain the 
following normal variational equation (NVE) 

V"{r) = -xHt) r]{T), xi(r) = kt , (2.4) 

which is solved by the parabolic cylinder functions I?-i/2((^ — 1)'*''") ^^'^ -^-i/2((^ + ^)kt). These 
are not Liouvillian functions and this implies the non-integrahility of the original Hamiltonian (2.3). 
To illustrate the previously discussed algebraization procedure let us consider the Hamiltonian 

H=\{pl+pl) + \{xl-xlf (2.5) 

which is equivalent to (2.3) by the trivial canonical transformation 

X1-X2 X1+X2 P1-P2 P1+P2 

x,^-^, .2^-^, Pi-^, P'^^TT- ^'-'^ 

Restricting to the invariant plane P = {{xi, X2,pi,P2) ■ X2 = P2 = 0} we obtain an integrable 
subsystem 

1 „ 

xi=Pi, pi = --x^. (2.7) 

The Hamiltonian vector field on this solution plane is X^^ = —^xf, X^^ = 0, X^^ = xi, X^^ = 0. 
The NVE is now obtained by expanding the Hamiltonian equations for the coordinates normal to 
the invariant plane, i.e. X2 and p2, along (2.7). Denoting 6x2 = rj we thus obtain 

f, = \x\{t)rj , -x,{t) = 23/4/.V4,n (2.8) 
where xi is a solution of (2.7) representing the curves in P parametrised by the Hamiltonian /i > 0. 
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We can algebrize this NVE by changing the variable t ^ x = xi (t). This leads to (' = ^) 

97.3 o™2 

rj" - ^ n' - , 1] = . (2.9) 

' 8h-x^' 8h-x^' ^ ^ 

Further changing the variable — > = {8h — x^^^^t] one obtains the NVE in the normal form 



^„ , 8hx'^ + 2x^ 
' {8h - x^y 



^"+ ... .4^2 ^ = 0- (2-10) 



This equation is solved in terms of hypergeometric functions, i.e. it has no Liouvillian solutions. 
This implies non-integrability of (2.5), in agreement with the previous discussion of of the equivalent 
system (2.3). 

To summarise, proving non-integrability of a Hamiltonian system can be achieved by demon- 
strating that given an integrable subsystem defined by some invariant plane in the phase space, the 
corresponding normal variational equations are not integrable in quadratures and thus do not admit 
sufficiently many conserved quantities. 



3 String motion in p-brane backgrounds 

In what follows we shall consider classical bosonic string motion in a curved background. In the 
conformal gauge the action is 

/ = y dadr [GMN{Y)daY^d''Y'' + e'''BMN{Y)daY''dbY''] , (3.1) 

and the corresponding equations of motion should be supplemented by the Virasoro constraints 

GmnY^'Y^ = 0, (3.2) 
Gmn (y^^y^ + Y^^Y^) = . (3.3) 

In this section we shall consider the case of Bmn = and the target space metric given by the 
following p-brane ansatz (Q > 0) 

ds2 = /-2(r) dx^'dx'^?]^^ + /2(r)(dr2 + r2 d^l) , (3.4) 
f{r)=(l + ^-Y\ n.m^Q (3.5) 



rj^y = diag(-l, +1, +1), = 0, 

where dil| is the metric on a /c-sphere. 
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For k = 8—p, n = 7—p > and m = j this metric describes the standard single p-brane geometry 
in = 10 supergravity. It may be supported by a R-R field strength and dilation backgrounds, but 
since they do not couple to the classical bosonic string we shall ignore them here. 

Keeping k, n, m,p generic also allows one to describe lower-dimensional backgrounds representing 
(up to a flat factor) some special brane intersection geometries. In particular, the metric (3.4) with 
nm = 1 describes a one-parameter interpolation [60, 61] between flat space (for Q = 0) and the 
AdSp+2 X S'^ space (for Q — ^ oo).^ 

Special cases include: 

{i) n = A, m = J, k = 5, p = 3: D3-brane interpolating between flat space and AdS^ x S^; 

(ii) n = 2, m = ^, k = 3, p = 1: D5-D1 (or NS5-F1 with non-zero Bmn, see next section) 
background [62,63] with Q5 = Qi = Q interpolating between flat lOd space and AdSs x x T^; 

(iii) n = 1, m = 1, k = 2, p = 0: four equal-charge D3-brane intersection [64] (or U-duality 
related backgrounds, see, e.g. [65]) interpolating between flat lOd space and AdS2 x S"^ x T^; this 
may be viewed as a generalised Bertotti-Robinson geometry. 

Below we shall show that the point-like string (i.e. geodesic) motion in this background is 
integrable. We shall then demonstrate that extended string motion in the p-brane geometry with 
p = 0, ..,6 and in the backgrounds (ii), (iii) is not integrable for generic values of Q, despite 
integrability being present in the limits Q = and Q = 00. 

3.1 Complete integrability of geodesic motion 

The symmetries of the metric (3.4) are shifts in the coordinates giving p+1 conserved quantities 
and spherical symmetry which gives k/2 conserved commuting angular momenta for even k and 
{k + l)/2 for odd k (generators of the Cartan subgroup of SO{k + 1)). Thus for k > 3 the spherical 
symmetry does not, a priori, provide us with sufficiently many conserved quantities. 

Parametrising the fc-sphere by /c + 1 embedding coordinates the effective Lagrangian for point-like 
string motion is 



The first integrals for and r are E^^ = f ^(r)x^ and the Hamiltonian respectively. The equation 
of motion for reads 



®More precisely, the metric (3.4) interpolates between the flat-space region (r" >• Q) and the AdSp+2 x S'^ region 
(r" <^ Q). Moreover, since f{r) = (c+ ;^)'" gives a solution for any constant c including c = 0, AdSp+2 X 5''= is also 
an exact solution. 



C = /-2(r)x^i;'^7?^, + f\r){r^ + r^y-^) + A{f - 1). 



(3.6) 



a,(/2(r)rV) + Ay' = 0. 



(3.7) 



7 



Taking the scalar product of this equation with y and with y and using = 1, = we 

conclude that 

A = f'^{r)r'^if = f~'^{r)r~'^w? , = f^{r)r^if = const . (3-8) 

Substituting this expression into (3.7) allows one to eliminate A giving 

drif {ryy') + m^f-\r)r-Y = . (3.9) 

Multiplying this equation by f'^{r)r'^y^ for fixed i and integrating once we obtain the constants of 
motion C* (no summation over i) 

f\r)rHy'? + m\f)^ = C\ ^C' = 2m' . (3.10) 

i 

We can rewrite the in terms of phase space coordinates (y*, = 2/^(r)r^y*) as 

C^ = l{p'f+m\y^f . (3.11) 

The integrals C* are easily checked to be in involution. Only k of the k + 1 integrals C* are 
functionally independent and correspond to the angles on the sphere. Altogether with the p+1 
constants for the coordinates and the Hamiltonian this gives in total k+p + 2 constants of motion 
which is the same as the number of degrees of freedom. 



3.2 Non-integrability of string motion 

Let us now study extended string motion in the geometry (3. 4), (3. 5). We shall choose a particular 
"pulsating string" ansatz for the dependence of the string coordinates on the world-sheet directions 
(r, cr): we shall assume that (i) only x", r and two angles (p, 6 ol S"^ C (with dJl^ = + 
s\v? 9 dcjP') are non-constant, and (ii) x^,r,d depend only on r while (f) depends only on a, i.e. 

xO = t(r), r = r(r), ^ = ^a), 6 = e{T) . (3.12) 



^Note that when only one of the angles of S'' is non-constant, 1-d truncations are not sufficient to establish 
non-integrability, see below. 



8 



This ansatz is consistent with the conformal-gauge string equations of motion and the Virasoro 
constraint (3.2). The remaining string equations of motion and the Virasoro constraint give 



i = Ef , E = const , (3.13) 

(j) = u = const , (j) = i/a , (3.14) 

2dr(,fr) = dr{f)E' + dr{f)r^ + dr{r^f){e^ - sin2 e) , (3.15) 

drifr^e) = - fr"^!^"^ sine cos , (3.16) 

E^ = r'^ +r'^e^ + u'^r^sm^e . (3.17) 

Thus the string is wrapped (with winding number v) on a circle of whose position in r and 6 

changes with time. The equations for r and 9 can be derived from the fohowing effective Lagrangian 

C = f (^r2 + r2^2 _ ^,2^2 g.j^2 0^E^') , (3.18) 



with the corresponding Hamiltonian restricted to be zero by (3.17). We shaU show that this Id 
Hamiltonian system is not integrable, implying non-integrability of string motion in the p-brane 
background (3.4). 

Let us choose as an invariant plane {{r,6;pr,pe) '■ ^ = fj P6» = 0}, corresponding to the a string 
wrapped on the equator of and moving only in r. Then (3.17) gives + i^^r^ = E'^ which is 
readily solved by (assuming e.g. that r(0) = 0) 

E 

r = r(T) = — sin(z^T) . (3.19) 

According to the general method described in section 2, we have to show that small fluctuations 
near this special solution are not integrable, or, more precisely, that the variational equation normal 
to this surface of solutions parametrized by E and has no Liouvillian solutions along the curves 
inside the surface. 

For the invariant subspace {6 = ^, pg = 0, r = f, pr = 2/2(f)f} the Hamiltonian vector field is 
= f, =Pr = ^ff'f"^ + 2/^f, = 0, XP" = (3.20) 

and thus the normal direction to this plane is along 9 and pg . Expanding the Hamiltonian equation 
for 9 around = f, ^ = 0, r = f one obtains the NVE (69 = tj) 

■■ , o \f-i , ^ ■ ^ x/'(f sin(i/T))\ . 2 n ron-,\ 

r? + 2i/cot(i^r)( 1 + -sm(z/T)— |— - — —jV-^ = . (3.21) 
V /(-sm(z^r))/ 

If /'// is a rational function (as is the case for (3.5)) this equation can be algebrized through the 
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change of variable r — )• x = sin(z^T) giving 



^ + ( 3 - + 2-7^^ K' - ^-f ^ = 0- (3-22) 



X 1 — 1/ /(-px) ^ 1 — x'^ 

Using the expUcit form of /(r) in (3.5) this NVE reduces to 

r] +{--- 0-2—. Ai]-- ^r? = 0. 3.23) 

' \x l-x2 x(i/"Q + S"x")/ l-x2' ^ ' 

We shah assume that v ^ Q a& otherwise the string is point-hke and we go back to the case of 
geodesic motion. 

Bringing (3.23) to the normal form by changing the variable to ^(x) = g{x)r]{x), where ^(x) is a 
suitably chosen function, one obtains 

r 2 + x^ (m — l)mn? mn(n — 1 — (n — 2)x^)i _ 

^ ^ U(x2 - 1)2 ~ x2(l + 5x")2 ^ x2(x2 - 1)(1 + Bx^) K = ^ ' ^ ' 

There are two special limits of (3.23) and (3.24) 

Q^O: ^" + ^(^^^^ = 0, (3.25) 

^„ 2x2 + x'^-4(mn)2(x2-l)2 + 4mn(l-3x2 + 2x'*)^ ^ ,^ 
g^oo: ^+ 4x2(x2 - 1)2 ^ = 0' (^-^e) 

corresponding to the flat space-time and AdS2 x 5*2 for nm = 1. 

We can now apply the Kovacic algorithm^ to determine whether these NVEs do not admit 
Liouvillian solutions and thus the identity component of their Galois group is not solvable which 
would imply non-integrability. 

For the special cases Q — t- and for Q — t- oo with nm = 1 Liouvillian solutions are found which 
is consistent with the fact that string motion in flat space-time and in AdS2 x 5*2 is integrable. 
However, for a finite value of Q Liouvillian solutions do not exist for generic values of E, v for the 
cases of a p-brane background with p = 0, ..,6 and the intersecting brane backgrounds of (ii) and 
(iii) mentioned above. This implies non-integrability of string motion in those special cases of the 
general background (3. 4), (3. 5). 

Finally, let us comment on the special case of a 7-brane in 10 dimensions when the transverse 
space is 2-dimensional, i.e. the string motion is described by 

C = r\r)dax^d''x^r,^, + fir) (S^raV + r^daOd'^e) . (3.27) 



We use Maple's function kovacicsols. 
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It is easy to see that if we truncate this model to a 1-parameter system by taking each of the 
string coordinates to be a function of r or o" only then the corresponding equations always admit 
constants of motion for and 9, and the Virasoro condition then provides the solution for r. 
To test integrability in this case one is to consider some more non-trivial truncations. Assuming 
the spatial coordinates are constant the corresponding effective metric is 3-dimensional: ds'^ = 
— f~'^{r)dt^ + f'^{r){dr'^ + r'^dO'^). Integrability of string motion in such a background deserves 
further study. 



4 String motion in NS5-F1 background 

Let us now include the possibility of a non-zero Bmn coupling and consider the case of string motion 
in a background produced by fundamental strings delocalised inside NS5 branes [3,62,66,67] (here 
we use the notation = x^ = z and i,j = 2, 5) 

ds^ = H^^{r){-dt^ + dz^) + dxMxi + H^{r){dr'^ + r^d^i) , (4.1) 
dO| = de"^ + sin^ e d(/j2 + cos^ 6 d<i? , 
B = -H{^{r) dtAdz + Q5 sin^ 9 dip A d(l) , (4.2) 

^5 = 1 + ^, Hi = l + ^, (4.3) 

with the dilation given by e~^* = . The corresponding string Lagrangian in (3.1) C = {Gmn + 
BMN)d+Y^^d-Y^ (where d± = dQ±di) is 

C = -T. — d+u d-v + d+x^d-Xi -\ 7, — d+rd-r 

+ {r"^ + Q5) {d+ed-9 + sin^ 6 d+^d-ip + cos^ 6 a+(/>9_(/>) + Q5 sin^ 9 {d+ipd-c/) - d-ipd+c/)) , (4.4) 
u = —t + Z , V = t + z . 

It interpolates between the flat space model for Qi,Q5 = and the SL{2) x SU{2) WZW model 
(plus 4 free directions) for Qi, Qs 00. Both of these limits are obviously integrable. 

Another integrable special case is Qi = 0, — >■ 00 when we get the SU{2) WZW model plus 
flat directions. 

The opposite limit of = 0, Qi — ?• 00 is described (after a rescaling of coordinates) by £ = 
9_i_M d-v + d^rd-r + r'^{d^9d-9 + sin^ 9 d^ipd-ip + cos^ 9 d^(t>d-4>)+ free directions. Solving the 
equations for u, v as (cr^ = r it o") 

r^d+u = f{a+) , r^d.v = f{a-) , (4.5) 
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one finds the effective Lagrangian for r being 



C = d+rd^r-^^^-^^^^^ + r^{d+ed^e + sm^ed+ipd-ip + cos'^ed+(l)d^(l)) . (4.6) 

As we shall find below, this model is not integrable.^ 

Let us note that in the limit — )• oo the action (4.4) reduces to a combination of the SU (2) 
WZW model, flat directions and the following 3-dimensional sigma model 

^ = 1 I ^ ^2o d+ud^v + Q^d+pd^p , p = Inr , (4.7) 

This model is T-dual to a pp-wave model with an exponential potential function [34,51] and is thus 
related to the SL{2) WZW model by a combination of T-dualities and a coordinate transformation 
(it can be interpreted, upon a rescaling of u, v, as an exactly marginal deformation of the SL{2, R) 
WZW model [35]). It should thus be integrable. Indeed, solving the equations for z it t as in (4.5) 
we get the following effective Lagrangian for p (cf. (4.6)) 

C = Q^d+pd.p - /(a-)/>+)(l + Qie-^O . (4.8) 

Since f{a~)f{a'^) can be made constant by conformal redefinitions of (reflecting residual gauge 
freedom in the conformal gauge) this model is equivalent to the Liouville theory, i.e. it is integrable. 
Indeed, the results of our analysis below are consistent with this conclusion. 



A point-like string does not couple to i?-field, so geodesic motion in the background (4.1) can 
be shown to be integrable in the same way as in the previous section. To study extended string 
motion let us consider the following ansatz describing the probe string being stretched along the 
fundamental string direction z (that may be assumed to be compactified to a circle) and the 
angle ip, i.e. 

t = t(r), z = z{a), x' = 0, r = r(r), = ^(r), ip = ip{a), 4> = 4>{r) . (4.9) 
Then the string equations are solved by 

i = KiiJi(r(T)) + K2 , z = K2a, Lp = va, (t>= or \u^r i \\ ' (^'^O) 

H(r)i?5(r(r)) 

^Let us note that the constant term in the harmonic function of the fundamental string background (1 in H\) can 
be changed by a combination of T-duality, coordinate shift and another T-duahty [51,61], so that the fundamental 
string background is not integrable for any value of Q\ if it is not integrable for large Q\. A similar shift in the 
harmonic function can also be achieved for D-brane backgrounds but this involves S-duality (mapping e.g. Fl to 
Dl) which is not a symmetry of the world-sheet string action; thus there is no contradiction with integrability of the 
D3-brane background at large Q. 
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The resulting Id subsystem of equations for r and 6 is described by the following effective Lagrangian 



C = H^{r) rVr^^^ - i/V^sin^^ - Qlu'^r^^H^^{r) cos'^ 9 + K?ifi(r) , 



(4.11) 



which should be supplemented by the consequence of the Virasoro constraint 

H5{r)\r'^ + + i^'^r'^ sin^ 9 + Qlu'^r-'^H^^ (r) cos"^ 9] - K;?i?i(r) = 2kiK2 , (4.12) 

implying that the Hamiltonian corresponding to (4.11) is equal to 2kiK2- This system admits the 
special solution 



2 



r = f{T) , H^{r){f^ + v^f^) = 2kiK2 + K{Hi{f) 



(4.13) 



which may be chosen as an invariant plane in the phase space. The corresponding Hamiltonian 
vector field is 



= f, = 2dr{fH^{r)), = 0, = . 



(4.14) 



Expanding along the normal direction to this solution plane we find for the NVE (^ = f + ??) 



,Hi{r) 2kiK2 



+ 



2-2 

v r 



Q 



Changing the variable r to a; = f(r) we obtain the NVE in the normal form (' = ^) 

fix) + U{x)ax) = , 



(4.15) 



(4.16) 



where 



=4^ 



Q5iQ5-2x^ 



+ 



4Q5a;4-2a;6 + 2Qig2(2Q5 + 2;2 



+ 



(Q5+x2)2 ^Q^+x^){-Qiq^ + x^[Q5+x^ -q^{l + a)]) 
3{Qiq^+xy 



iQiq^+x^[-Q5-x'' + q''{l + a)]f 







q = — 

V 



a^2^. (4.17) 



One finds that already for a = 0, i.e. for K2 = 0, this NVE has no Liouvillian solutions for general 
values of Qi and Q^, implying non-integrability of string motion in the NS5-F1 background with 
generic values of the charges. Our results are summarised in the table below. 

Let us discuss some special cases. Taking Qi = or Qs = one can absorb the remaining brane 
charge by rescaling x, q and k\. The resulting equations do not admit Liouvillian solutions for 
arbitrary values of q which implies non-integrability of string motion in the NS5-brane or in the 
fundamental string backgrounds. 
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In the limit Qi — )• oo, — t- oo, in which the non-trivial part of the string action becomes that of 
the SL(2, R) x SU{2) WZW model, one finds Liouvillian solutions in agreement with the expected 
integrability. The same applies to the case of Qi — s- 0, — )• (X) described by the SU(2) WZW 
model plus free fields and, in fact, to the model with — )• oo and arbitrary Qi described by 
(4.7)(4.8). 

In the opposite case of — >• 0, Qi — )• oo described by (4.6) integrability appears to be absent as 
we did not find Liouvillian solutions (see table). ^'^ 

Since 2d duality transformations of string coordinates preserve (non) integrability, similar conclu- 
sions can be reached for string backgrounds related to this NS5-F1 background or to its limits via 
T-dualities (and coordinate transformations). In particular, this applies to the pp-wave background 
related to the fundamental string by T-duality [51] which is studied in in Appendix B. 



5 Concluding remarks 

We have shown that for various p-brane backgrounds, for which the string sigma model interpo- 
lates between integrable flat and coset or WZW models and which has integrable geodesies, the 
corresponding extended classical string motion is not integrable in general. 

To demonstrate non-integrability we considered particular extended string motion for which the 
dynamics reduces to an effective one-dimensional system. The latter has special integrable subclasses 
of solutions, but perturbing near them leads to linear second-order differential equations that cannot 
be solved in quadratures. 

It would be interesting to understand better why switching on a non-zero D3-brane charge or 
moving away from the "throat" {AdS^ x S^) region leads to a breakdown of string integrability. 

Together with similar previous results [6,7,49] this supports the expectation that integrability of 
classical string motion is a rare phenomenon. String integrability is thus a much more restrictive 
constraint than integrability of particle motion. Still, in the absence of a general classification of 
integrable 2d sigma models one cannot rule out the possibility that there are still interesting exam- 
ples of integrable backgrounds (that cannot be obtained from known solvable cases by coordinate 
transformations combined with T-dualities, cf. [29,30]), that remain to be discovered. 



^ Let us note that the results about the existence of Liouvilhan solutions turn out to be independent of whether the 
truncated solutions contains any contribution from the B-flux of the fundamental string solution (a 7^ or a = 0), 
reflecting a special choice of our ansatz (4.9). 
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limit 



U{x) in NVE + U{x)C{x) =0, 7 = §7 



rescaling 



Qi ^ 

00 Qi — 00 

00 Qi — 

Qi — )■ 00 

00 Qi / 



Qs / Qi ^ 00 

Qs ^ Qi / 
Qs / Qi ^ 



4(g^(l+a)-a;^)^ 

(7-g^a)(73:-^+g^(2-3;^a) 
4(7a;^-g^(l+a;^a))^ 

1 

(7-g^(l+«))(7-r^-g^(-2+a:^(l+a))) 
4(7X^-g^(l+a;^(l+a)))2 

7-2g*(a;2 - 2) + 27-19^(1 + x'^f{l + Sx^) - 27-lg''(l + Sx^ + x'')a 
+ x\(l + x^f - + x^f + q\2x^ - l)a)] 

X [4(1 + x^f[x^ +x^ - <J^{7"^ + x^a)f'^ 



Ot — 7- QlOi 



x^+g-*(l-2a:^(l+a))+29^x^(5+x^(l+a)) 
4a:^(a;''-g^(l+3;^(l+a)))^ 

(l+a^)'t+2g^(x^-l)(l+x^)^(l+a)-q4(2x^-l)(l+a)^ 
4x2(l+a;2)2(l+x2-g2(i+a))2 



a; 

Kl,2 

a; 

Kl,2 
Kl,2 

fi;i,2 



/QlKl,2 



/Q5'tl,2 
/QlKl,2 



\/^Kl,2 



limit 



spacetime 



Ri''^ X 53 



Fl 

NS5 



2-form contributes to the 
truncated system 



a = 



a / 



NVE has Liouvillian 
solutions 



a = 



a 7^ 







Q5 

Qs ^0 



Qi ^ 
00 Qi — )■ 00 
00 Qi 
Qi — 7- 00 
00 Qi / 
Qi 00 
Qi 

Qi ^0 



X 
/ 
/ 
X 

/ 
/ 

X 

/ 



X 

/ 
/ 
/ 
/ 
/ 
/ 
/ 



/ 
/ 
/ 

X 

/ 

X 
X 
X 



/ 
/ 
/ 

X 

/ 

X 
X 
X 
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A Integrability and differential Galois theory 

Integrability in the sense of Liouville implies that the equations of motion are solvable in quadra- 
tures, i.e. the solutions are combinations of integrals of rational functions, exponentials, logarithms 
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and algebraic expressions. Solutions of this type are referred to as Liouvillian functions. Thus the 
question of integrability is related to the question of when differential equations admit only Liou- 
villian solutions. This is answered by differential Galois theory and the main integrability theorem, 
which our analysis relies on, can be stated as follows [69]: 

Theorem Let H be a Hamiltonian defined on a phase space manifold M with an associated 
Hamiltonian vector field X. Let P he a submanifold of M which is invariant under the flow of 
X with the invariant curves T and let X\p be completely integrable. Also let G be the differential 
Galois group of the variational equation of the flow of X\p along T normal to P. If X is completely 
integrable on M then the largest connected algebraic subgroup of G which contains the identity is 
abelian. 

Since differential Galois theory is not frequently used in physics we shall give a brief introduction 
and illustrate some basic concepts on a simple example below. We shall follow ref. [70]. First, let 
us give some basic definitions: 

1. A differential field is a field F with a derivation Dp, which is an additive map Dp : F ^ F 
thus satisfying Dp(ab) = DF{a)b + aDp[b) , Va, b & F. 

2. A differential homomorphism/isomorphism between two differential fields Fi and F2 is a ho- 
momorphism/isomorphism / : Fi — )• i*2 which satisfies Dp^[f{a)) = f{Dp^{a)) , Va G Fi. 

3. A linear differential operator L is defined by 

L{y) = + a„_i?y("-i) + ... + aoy (A.l) 

where y*^") = Dp{y) and Oj G F. 

4. A differential field E with derivation Dp is called a differential field extension of a differential 
field F with derivation Dp iS E ^ F and the restriction of Dp to F coincides with Dp. 

5. Elements a £ F whose derivative vanishes are called constants of F. The subfield of constants 
of F is denoted by Ci;'. 

In order to investigate the types of solutions a differential equation admits one has to encode the 
relations and symmetries of the independent solutions. Relations can be encoded in the Picard- 
Vessiot extension of a differential field which includes the solutions of L{y) = 0: 

Definition An extension E ^ F is called a Picard-Vessiot extension of F for L{y) = iff: 
E is generated over F as a differential field by the set of solutions of L{y) = in 
E and F have the same constants; 

L{y) = has n solutions in E linearly independent over the constants. 
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One can show that the Picard-Vessiot extension always exists for differential fields with an alge- 
braically closed field of constants and this extension is unique up to isomorphisms. 

Symmetries of the solution space are linear transformations which, applied to a solution, give a 
new solution. The simplest formulation is in terms of automorphisms of the Picard-Vessiot extension 
and leads to the notion of the differential Galois group: 

Definition For a differential field F with the Picard-Vessiot extension E ^ F the differential Galois 
group G{E/F) is the group of differential automorphisms a : E ^ E whose restriction to F is the 
identity map, i.e. 

G{E/F) ={a: E^ E\ a{a) = a Va e F} . 

A solution y{x) of a linear differential equation L{y) = is called 
algebraic if it is a root of a polynomial over F; 
primitive if y'{x) G F, i.e. y{x) = f{z) dz for f{x) G F; 
exponential if € F, i.e. y{x) = exp ( J'^ f{z) dz) with f{z) € F. 

Definition The solution y{x) is called Liouvillian if there exist differential extensions Ei, i = 0, ..,m 

F = Eo C El C ... C Efn = E 

such that y{x) = ym{x) £ E, Ei = and yi{x) is algebraic, primitive or exponential over 

E-i-i- 

We will be interested in second order homogeneous linear ordinary differential equations 

y"{x) + a{x)y'{x) + b{x)y{x) = , (A.2) 

where a{x),b{x) € F and F = C{x) is the differential field of rational functions over complex 
numbers. The equation (A.2) can be cast into the normal form 

y"{x) + U{x)y{x)=0 , (A.3) 
U{x) = -^a^{x) - ^a'{x) + b{x) , U{x) G F (A.4) 

by the transformation y{x) — )■ y{x) exp ( — ^ a{z)dz^ . 

Let us identify the corresponding Galois group. Let E he a Picard-Vessiot extension oi F = C{x) 
and zi{x), Z2{x) £ E he two independent solutions of (A.3). The Galois group G{E/F) consists 
of differential automorphisms a which act on the space of solutions of (A.3), i.e. map solutions to 
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solutions. Thus 



a{zi{x)) = CijZj{x) , i,j = 1,2 . (A.5) 

I we obtain 

z[{x) / 

a{w{zi{x), Z2{x))) = detC w{zi{x), Z2{x)) (A. 6) 

where detC is the determinant of the Cij matrix in (A.5). The Wronskian is a constant since 
w'[x) = Zi{x)z2{x) — Z2{x)z'l{x) = by virtue of (A. 3), i.e. w{x) G F. Therefore, (t{w{x)) = w{x), 
giving detC = 1. Thus the Galois group for (A. 3) is a subgroup of SL(2,C). For more general 
equations the Galois group is a subgroup of GL(77,,C). 

Let us now address the main question: whether a differential equation can be solved in terms 
of Liouvillian functions. Since the minimal differential extension generated by all solutions is the 
Picard-Vessiot extension we therefore require that the Picard-Vessiot extension is generated by a 
tower of Liouvillian extensions. By the fundamental theorem of differential Galois theory this tower 
of extensions is related to subgroups of the Galois group through a bijective correspondence. One 
can show that for the Picard-Vessiot extension to be Liouvillian, the identity- connected component 
of the Galois group has to be solvable (see theorem 25 in [56]).^^ 

A simple example is provided by the equation y"{x) + 2xy'{x) = which has the normal form 

y"{x) - {x^ + l)y{x) = 0. (A.7) 



For this equation F = C(x) is the field of rational functions over complex numbers. The solution 
space may be represented as C(e2^^) © C(e2^^ J e~^^). The Picard-Vessiot extension is E' = 
Ox, 6 2^^, e~^^, / e~^'^) and the tower of extensions of E is 



Then the Galois group G{E/F) and its identity component are 

G = G0 = {( " ° ), a/0,c Gcj (A.8) 
This is a solvable group since it is a subgroup of SL{2, C) whose algebraic subgroups are all solvable 



^^A group is solvable if it has a subnormal series whose factor groups are all abelian, that is, if there are subgroups 
{1} < Gi < ... < Gn ~ G such that Gk-i is normal in Gk (i.e. invariant under conjugation) and Gk/Gk-i is an 
abelian group for all k = 1, n. 
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except for SL{2, C) itself. 



B String motion in pp-wave background 

Here we shall study string motion in a pp-wave metric (i = 1, ..,d) 

ds'^ = dudv + H{u,x)du'^ + dx'dxi . (B.l) 

In conformal gauge the equation for v reads d^^d-U = and is solved by n = /(fj"*") + /(cr~). We can 
fix the residual conformal symmetry by choosing u = pr. Then v is determined from the Virasoro 
constraints and we get the following effective Lagrangian for Xi 

£ = x'x' - x'£ + p'^H{u, x), (B.2) 

which describes light-cone gauge string motion in a potential. 

One familiar example is the Ricci flat space with H{x) = ^ijX^x^ , ^\ = 0. Another is the pp-wave 
limit of the AdS^ x background for which H[x) = x^Xj [44]. Here the string motion takes place 
in a quadratic potential and is obviously integrable. Penrose limits of the brane backgrounds we 
consider in this paper also take the pp-wave metric form with H{u, x) = hij{u)x^x^ (see e.g. [46]). In 
the light-cone gauge string motion takes place in a quadratic potential with a r-dependent coefficient 
and it is solvable [47] but formally the resulting model is not integrable. 

There are also other integrable examples with H(x) corresponding, for example, to the Liouville 
or Toda potential. 

The case that we shall study below is the pp-wave solution T-dual to the fundamental string 
]51,66] for which H is a. harmonic function 

H{x) = 1 + ^ , = x'xi . (B.3) 

Not surprisingly, the conclusions will be the same as for the fundamental string in section 4: the 
geodesic motion is integrable but an extended string motion is not. 

In the point-like string limit (B.2) takes the form (we set rr* = ry* with = 1) 

£ = r2 + r2^2^p2^(r)+A(y2_l) . (B.4) 

The corresponding geodesic motion is completely integrable since we have d—1 constants of motion 
from the coordinates (see section 3.1) plus the Hamiltonian. 

^■^For the relevant case oi h ^ ^ the Kovacic algorithm gives no Liouvillian solutions for string modes that depend 
on (T. 
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Let us now consider an extended string moving only in a 3-space dx^dxi = (ir^+r^((i0^+sin^ 9 dcf)^) 
and choose the foUowing ansatz (similar to the one in (3.12)) 

r = r{T), (j) = (l){a) = ua, 6 = 0(r) . (B.5) 

Then the effective Lagrangian for r and 9 takes the form 

£ = r2 + _ sin^ 9 + p^H{r) . (B.6) 

If we also assume that v = (i.e. the string also moves in the light-cone direction) then the Virasoro 
condition for (B.l) is equivalent to the vanishing of the Hamiltonian corresponding to (B.6) 

^2 _^ ^2^2 _^ ^2^2 g-j^2 g _ p2fj^^^^ ^ . (B.7) 

Restricting motion to the invariant plane {{r,pr,9,po) : 9 = ^, pg = 0} gives a one-dimensional 
integrable system parametrised by v and p with the Hamiltonian as the constant of motion. Imposing 
the zero Hamiltonian constraint (B.7) gives the solution for r 

r2 + j,2^2 ^ p-2H{r) , r = f{T) . (B.8) 

The Hamiltonian vector field on the plane of these solutions is 

X'' = f, X'P- =2f = -2i/V + p^H'{f), = 0, = (B.9) 

and the direction normal to this plane is along 9 and pg. Expanding the equation of motion for 9 
gives the NVE {r^ = 59 = 9 - ^) 

ri + 2^7) - u^rj = . (B.IO) 
r 

Changing the independent variable t ^ r = f{T) one obtains (' = ^) 

„ /2 \p^H'{r) - z/Vn , 1/2 



V + -+ 2m\ 22 h-^T7r^ — r^^ = o 

\r p^H(r) — i/^r^ / p'^H[r) — z/^r^ 
Bringing this equation to the normal form via the change of variable r/(r) = g{r)^{r) we get 

e + U{r)i = Q, (B.12) 

^ _-{d- 2)q^Q [{d - 6)Qr^ + 4{d - 3)r'^] + Aq^r'^ [{d^ -2d + 2)Qr^ + 2r'^] + 4r2+2'^ _ ^ 

16[g2(Q^2 + ^d) _^2+d]2 ~ ' 
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where q=^- For 2 < d < 13 the Kovacic algorithm does not yield LiouvilHan solutions for general 
finite non-zero values of Q and q which implies non-integrability. 

Note that in the limit Q — )• 0, when the metric becomes flat, we get 

which has Liouvillian solutions in agreement with with flat space integrability. 
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